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ABSTRACT 

This paper investigates numerically influence of linewidth enhancement factor and injection current on the 
type of the route-to-chaos and associated operation states of laser diodes under external optical feedback. The study 
is based on numerical solutions of a time-delay model of rate equations, and the solutions are employed to construct 
bifurcation diagrams and to examine the output time fluctuations and phase portraits of the laser. The simulation 
results show that linewidth enhancement factor causes significant changes in the route-to-chaos and the laser states. 
The state of the laser is identified into six distinct regimes, namely, continuous wave, periodic oscillation, period 
doubling or two, period-three, period-four oscillations and chaos, which is depending on the value of the linewidth 
enhancement factor, feedback strength and injection current level. The route-to-chaos is period-four when linewidth 
enhancement factor is less than or equal to three. The route is period-three when linewidth enhancement factor 
increases up to 7. The feedback strength when the laser transits from continuous wave to periodic oscillation, period 
doubling or chaos state decreases with the increase in the linewidth enhancement factor. Decreasing the injection 
current and increasing the linewidth enhancement factor stimulating the laser to operate in continuous and periodic 
oscillations. 
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1. INTRODUCTION 

Laser diodes (LDs) such as InGaAsP/InP lasers emitting at 1.2- 1.6 pm are widely used as light sources 
for long-distance optical communication systems [1]. The operation states and route-to-chaos of such lasers are 
strongly modified by optical feedback (OFB) from an external reflector. The variation of the operation states and 
types of the route-to-chaos of LDs depend on various factors, which divided into (i) internal parameters, such as 
nonlinear gain (NLG) and nonradiative recombination life time (NRRLT); and others, and (ii) external parameters, 
such as OFB strength, external cavity length and injection current. While weak OFB levels result in linewidth 
narrowing and improved frequency stability, moderate to strong OFB levels might induce the laser to switch to a 
state of significant spectral broadening and dynamical complexity, which has been called “coherence collapsed 
state” or chaos [1-4]. 

However, the operation and states characteristics of the LD are influenced by the linewidth enhancement 
factor (LEF), also known as a-factor, which has a great importance for LDs, as it is one of the main features that 
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distinguish the behavior of LDs with respect to other types of lasers. The LEF influences several fundamental aspects of 
LDs, such as the linewidth, the chirp under current modulation, the mode stability, the occurrence of filamentation in 
broad-area devices. In synthesis, the dynamics of LDs is greatly influenced by the LEF, which is of particular interest for 
the study of injection phenomena, OFB effects, and mode coupling as occurring in VCSELs [5-14]. 

The route-to-chaos is one of the interesting nonlinear phenomena associated with OFB; it may include dynamics 
with period-doubling (PD), sub-harmonics (SH) or quasi-periodicity (QP) [15-17]. PD is observed in the limit of short- 
external cavity [16, 17]; it follows a region of periodic oscillation (PO) and is characterized by the PO frequency and its 
half harmonic. Kao et al. [16] predicted that this frequency should correspond to Hopf-bifurcation (HB) point (the onset of 
PO). SH transitions to chaos were predicted as an intermediate state from PD to QP associated with increasing the external- 
cavity length [16]. It is characterized by the HB frequency and one of its rational SH, and it is followed by mode locking 
with these frequencies [16]. Ahmed et al. [18] used the ratio of the relaxation frequency to the external-cavity resonance 
frequency to classify the route-to-chaos in LDs. They showed that the route is PD when the frequency ratio is less than 
unity. The route becomes SH with the periodic oscillation frequency when the ratio is higher than unity and less than 2.25. 
When the ratio exceeds 2.3, the route is QP characterized by the compound-cavity frequency and the relaxation frequency 
as well as their frequency difference [18]. Abdulrhmann et al. [19-21] investigated the influence of NLG and NRRLT on 
the operation states, route-to-chaos and noise of LDs. He showed that, the value of the OFB rate at which the transition 
from CW to PO or chaos state occur increases, as the intensity of the NLG/NRRLT is increased/decreased. Under strong 
OFB, the increase/decrease in the NLG/NRRLT causes change of chaotic dynamics to PO or CW operation depending on 
OFB strength [19, 20]. The relative intensity noise is found to be as low as the quantum noise level when the laser is 
operated in CW or PO region and at relatively large/small values of the NLG/NRRLT [21]. 

However, based on our knowledge, there is no report that studied the effect of LEF, OFB strength and injection 
current ratio on the route-to-chaos characteristics and associated states of operation for LDs. In this article, we apply the 
time -delay rate equation model in [18, 22] of OFB in LDs to study influence of LEF on the route-to-chaos and associated 
laser operation. In such a model, OFB is treated as time delay of the laser radiation due to multiple roundtrips between the 
laser front facet and the external mirror. Intensive computer simulations are run to investigate the operation and dynamics 
of 1550-nm InGaAsP/InP LD with a short external cavity using a wide range of injection current. The route-to-chaos and 
associated laser operation is classified in terms of the bifurcation diagrams of the peak value of the photon number at each 
feedback strength and injection current. We identified the operation states of LD into six distinct operating regimes, 
namely, continuous wave (CW), PO, PD, period-three (3xP) and period-four (4xP) route-to-chaos and coherence collapse 
or chaos depending on the value of the LEF, OFB strength and injection current ratio. The results show that, variation of 
LEF cause significant changes in the route-to-chaos and operation states of the laser. The value of the OFB strength at 
which the transition from CW to PO or chaos state occur decreases, as the value of the LEF is increased. Under lower 
values of injection current and by increasing LEF, the laser operation states are only CW and PO states and there is no 
chaos state. At lower values of the injection current, highest LEF works in such a way to induce CW and PO operations. 

The paper is organized as follows; in sec. 2 the time-delay rate equation model of LD with OFB is introduced. 
These rate equations are solved numerically in Sec. 3. In Sec. 4, the bifurcation diagrams of the photon number of LD over 
wide ranges of LEF and injection current are shown. The time variations of the photon numbers and phase portrait of the 
induced route-to-chaos states of operation are also presented. The conclusions are presented in Sec. 5. 
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In this paper, we apply the model in [18] to investigate the effect of the LEF on the operations, and route-to-chaos 
of InGaAsP/InP laser over wide ranges of the OFB strength and injection current. The time-delay rate equations of the 
photon number S( t), optical phase (f t) and carrier number N( t), which describing the LD dynamics are formulated as 
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where ag(N-N,,)/V is the linear gain coefficient with a and N„ as material constants and g as the confinement factor of the 
optical field into the active region of volume V. a is the linewidth enhancement factor (LEF), and I is the injection current 
and N is the corresponding time averaged electron number. G,ho is the threshold gain level, which is determined by the 
loss coefficient k of the laser and mirror loss: 
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Various models have been advanced to account for gain nonlinearities based on the density-matrix theory [24-28]. 
In this paper, the nonlinear phenomenon of NLG is included in the rate equations in terms of coefficient B , which is given 
in the third-order perturbation theory of gain by [23, 24]. 
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where R IV is the dipole moment, T m is the intra-band relaxation time and N s is an injection level characterizing the NLG. 
The electron lifetime which should account for the non-radiative recombination by Auger processes since they have an 
effective contribution to the lasing mechanism in long wavelength laser is T s , which is defined as [25] 

- = A nr + B r (N/V) + C AUG (N/V) 2 , (6) 


where A m , is a coefficient for the nonradiative contribution arising from crystal impurities, B r is the radiative decay 
coefficient, and C A uo is a decay coefficient due to Auger processes [25]. The complex coefficient T describes the influence 
of OFB on the threshold conditions are given by [18, 22]: 
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where, m is an index for the round-trip, ip is a phase term combining the phase changes due to reflection by the external 
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mirror, R ex , n ex and L ex are the external mirror power reflectivity, refractive index and length of the external cavity, 
respectively. The combined phase y/ is defined in terms of the phase retarded due to reflection by the external mirror <j) ex , 
front facet, </> f , and due to a round-trip in the fiber cavity cot, as y/= <j) ex + <j) f + COT , where 0) is the emission circular frequency 
and r = 2 n ex L e Jc is the round-trip time. The strength of OFB is measured by the coefficient K ex , which is determined by 

the ratio of the external reflectivity R ex to the front -facet reflectivity R f as. 


= d’/Jr 


( 8 ) 


where 77 is the coupling ratio of the injected light into the laser cavity. The argument <j) of the complex OFB function T is 
obviously given by: 

^ = -tan~ 1 {lm[r]/Re[7’]] + /;r (9) 


where £ is an integer. Determining the value of <p in the two-dimensional space (Re[T]-Im[TJ) depends on both the signs 
and magnitudes of Re[T] and Im[T], 

3. PROCEDURES OF NUMERICAL CALCULATIONS 

Equations (l)-(3) are solved numerically with the fourth-order Runge-Kutta method. The simulation model is 
applied to InGaAsP/InP lasers emitting in the wavelength of 1550 nrn. The time step of integration is set as At =5 ps, which 
is short enough to get good resolution of the time variation of the photon number. The integration was carried out over a 
period of 6 ps, which is long enough to achieve stable dynamics of the laser. The injection current is varied between 1.05 
I,i l0 and 2.0 I th0 , where I th0 = 10.8 niA is the threshold current of the solitary laser. The influence of OFB is taken in terms of 
the strength of OFB K ex . The integration was first made without OFB (the case of the solitary laser) from time t = 0 until 
the first round trip time T, i.e., putting K ex — 0 and (p = 0 in the rate equations. The calculated values of S(t) and 0(t) are 
then stored for use as time delayed values S(t-T) and d(t-x) for the further integration of rate equations (l)-(3) including the 
feedback terms. We include three round trips (m=5). The calculated values of S(t= 0 to T) and &, /=() to T ) are then stored 
for use as time delayed values S(t- t) and 6(t- t) for integration of the rate equations over the period ;=2 rto 3 r including 
OFB terms. Then the calculated values S(t-t), S(t-2f), 6(t-f), and 6(t-2 T) are used as time delayed values for integration 
over the period t= 2 T to 5 T. Starting from this time, the integration is proceeded over a long period of time 7=4-6 ps by 
considering all terms of S(t-mT) and Oit-m T), m=0,l,2, 3, as time delayed values. The value of the combined phase change 
tj/ is set to be y/= 0.0, and the value of the external cavity length is set to be L ex — 0.03 m in the present calculations. The 
phase <p of the feedback light was chosen to vary continuously for time evolution because the solution of arc tangent is 
limited in the range of -n!2 to /z/2 in the computer work. The averaged values Tp are set as zero in this paper. The applied 
numerical values of InGaAsP/Inp lasers are given in table 1 in Ref. [19]. 
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In order to understand the effect of the LEF on the behavior of the external cavity LDs, we simulate the laser 
dynamics over wide ranges of the LEF, the OFB strength (in terms of K ex ) and injection current ratio. 

4.1. Bifurcation Diagrams Analysis 

Investigation of the influence of LEF on the route-to-chaos of the laser is done by simulating the bifurcation 
diagrams of the photon numbers S(t) in terms of the OFB strength K ex and injection current ratio I/I,ho- The bifurcation 
diagrams, which constructed by plotting the peak(s) of the temporal trajectory of the photon number S(t) at each level of 
OFB strength K ex , and at three values of the injection current ratio, l/l,i,o = 1-05 (near the threshold), 1.5 (well above the 
threshold) and 2.0 (far from the threshold) and at four values of LEF, a= 1.0, 3.0, 5.0 and 7.0, are plotted in figures 1-3, 
respectively. The results are normalized by the time ensemble average S ■ In these figures, we are interested in exploring 
an overview of the LD operation over a wide range of the LEF, injection current and OFB strength. 

General features of these diagrams are as follows. Under very weak OFB, the solution of the rate equations is still 
stationary and the laser operates in CW operation for which no points are plotted in the figure. With the increase in OFB 
strength, K ex , this stationary solution bifurcates first into a stable limiting cycle characterizing periodic oscillation (un- 
damped relaxation oscillation), which is represented by a single point in the diagram. The starting point of the bifurcation 
is called a Hopf-bifurcation (HB) point and is characterized by frequency f p , which is not necessarily equal to the 
relaxation oscillation frequency/^, and as indicated in the Figures 1 (a to d) by first arrow (black arrow) is moved toward 
decreasing OFB strength K ex from 0.028 to 0.0064 by increasing LEF from 1.0 to 7.0, respectively. 

With further increase in OFB, the solution of the rate equations bifurcates into a torus (multiple bifurcation points) 
followed by a chaotic state, and then the laser is attracted again to CW operation. As shown in the figures 1-3, the torus is 
PD in which periodic oscillation bifurcates first into two branches, where the trajectory of S(t) has two peaks of different 
heights in every two successive periods, which indicated in the Figures 1 (a to d) by second arrow (blue arrow) in which 
PO bifurcates first into two branches. With the increase in K ex , the oscillation period is multiplied to more than twice and 
the laser is attracted to transition to coherence collapse or chaos, which indicated in the Figures 1 (a to d) by fourth arrow 
(red arrow). When LEF a= 1 and 3, the route-to-chaos is 4xP as indicated in the Figures 1 (a, and b) by third arrow (green 
arrow), respectively. In Figures l(c and d), and when LEF, a = 5 and 7, the route-to-chaos is 3xP (green arrow). Then by 
increasing the LEF the route-to-chaos of the laser changes from 4xP route-to-chaos to 3xP route-to-chaos. 

Effect of the injection current on the states and route-to-chaos of the LD at several values of LEF is shown in 
Figures 2 and 3. Well above the threshold current I th0 , at I/I, h0 = 1.5, the states and route-to-chaos is changed from 4xP to 
PO route-to-chaos when a= 1.0 and from 3xP to 4xP route-to-chaos when a= 7.0 as shown if Figures 2 (a, and d), 
respectively. By decreasing the injection current to I/I th0 = 1.05 (near threshold), the laser operates in CW and PO and there 
is no chaos state, especially at higher values of LEF as shown in Figures 3(a-d). With the decrease in I/I t ho, the laser 
operation changes to CW and PO at all values of the LEF, and then the laser is attracted to transition to CW operation by 
increasing the OFB strength. 

In the following subsection, we characterize the three types of route-to-chaos when I/I th0 = 2.0 in terms of the time 
variation of photon numbers S(t) and corresponding phase portraits of the state of the LD. 
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Figure 1: Bifurcation diagrams of a LD under OFB when I/I, h0 = 2.0 (a) a= 1.0, (b) a= 3.0, (c) a= 5.0 and (d) a= 7.0 






Figure 2: Bifurcation diagrams of a LD under OFB when I/I„ l0 = 1.5 (a) a= 1.0, (b) a= 3.0, (c) a= 5.0 and (d) a= 7.0 


Impact Factor (JCC): 6.2879 


NAAS Rating: 2.40 


Route-to-Chaos and Operations of Laser Diode under 
Optical Feedback with Linewidth Enhancement Factor 


25 



Figure 3: Bifurcation diagrams of a LD under OFB when I/I th0 = 1.05 (a) a- 1.0, (b) a= 3.0, (c) a= 5.0 and (d) a= 7.0 



Linewidth Enhancement Factor a 

Figure 4: Dependencies of the OFB strength K ex at HB ppoint on the 
LEF a at iinjection ccurrent rratio I/I th0 = 1.05, 1.5, and 2.0 

The feedback level when a PO appears at HB point calculated from the numerical simulation of the rate equations 
using bifurcation diagrams shown in figures 1-3, as a function of the LEF is shown in figure 4. Notice that they are almost 
nonlinear functions (exponential decay) of LEF, which means that increasing the LEF decreases the feedback level 
required to destabilize a PO, PD, ... or a chaos state. Figures 1-3 and 4, demonstrate that, decreasing the LEF renders a 
steady state more stable by increasing the feedback level above. By decreasing the LEF, the level of the OFB strength 
corresponding to a HB point and route-to-chaos increases, which leads to increasing in the laser threshold current. We 
believe that, from the result shown in Figures 1-3 and 4, the instability of the LD can be reduced by decreasing the LEF 
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and decreasing the injection current ratio to values near from the threshold and increasing OFB strength, which helps to 
design a laser diode with high static and dynamic performance. So, it is clear that from Figures 1-3 and 4, any increase in 
the LEF leads to increase in the threshold current [8], And any increase in the OFB strength reduces the total cavity loss of 
the laser diode [26], which caused to increase the photon lifetime and increase the threshold carrier density, which increase 
the laser threshold current. Thus, from the physical and practical point of view, any increase in the laser threshold current 
leads to increase in the laser stability. These results support our conclusion that, at higher values of LEF and when the 
injection current be near the threshold the laser stability is improved and moves toward stable operations. 

In the following sections, we characterize the operation states of the laser and the two types of route-to-chaos in 
terms of the time variation of photon number 5(f ), and corresponding phase portraits for various states of the laser. 

4.2. Time Variations of S(T) and Phase Portraits at Several States of the LD 

In this section, we examine the time variation of photon number from back facet 5( f) and phase portrait of LD of 
the operation states of LD till reach to chaos shown in Figures. 5 and 6 at I/I th0 =2.0 and at three values of the LEF a = 3.0, 
5.0 and 7.0. The examination is done over a wide range of the operation states (PO, PD, and 3xP, 4xP route-to-chaos, and 
chaos). 

The simulated time variations of the photon number 5(f), of PO operation of the laser at K ex = 0.02, 0.015, and 0.01 
and when a= 3.0, 5.0 and 7.0 are shown in Figures 5 (a, e, and i ). The time variation of the PO laser operation exhibits 
uniform fluctuation with period-one. Theses PO states are confirmed by plotting the corresponding phase portrait in 
Figures 6 (a, e, and i ), The phase portraits are characterized by limit cycle attractor (period-one oscillation). By increasing 
K ex to 0.027, 0.02 and 0.015, which represent the regions of the PD operation, the laser operation becomes more 
complicated and the time variation becomes with period-two, which could be considered a kind of PD oscillations as 
shown in Figures. 5 (b, f, and j). Figures 6 (b, f, and j) show the phase portrait of PD operations, which are characterized by 
limit cycle attractor (period-two oscillation). When K ex = 0.0284, 0.023, and 0.017, the system nonlinearity increases 
irregularities of 5(f) attracting the laser oscillation into route-to-chaos as shown in Figures 5(c, g, and k). As shown in 
Figures. 5 (c) the time fluctuations of 5(f), induce four peaks, which represent 4xP route to chaos and this is confirmed by 
the phase portrait shown in Figure. 6 (c), which is characterized by limit cycle (period-four) attractor. By increasing LEF to 
5.0 and 7.0 the route to chaos is changed to 3xP as shown in Figures 5 (g and k), the time fluctuations induce three peaks 
and the phase portrait is characterized by limit cycle (period-three) attractor as shown in Figures 6 (g and k). When 
(K ex = 0.03, 0.025, and 0.018), the laser state becomes more complicated and enter to coherence collapse or chaos regime. 
Figures 5 (d, h, and 1), shows non-uniform or random fluctuations of the time variation of 5(f), which is appeared as a 
chaotic attractor in the phase portrait shown in Figures 6 (d, h, and 1). 

CONCLUSIONS 

We investigated the influence of LEF on the states and route-to-chaos of LDs. The simulations were applied to 
1550-nm In GaAsP/InP lasers over a wide range of the injection current. The route-to-chaos and the associated laser 
operations are classified by the bifurcation diagrams of the photon number, time-variations of the photon number. The 
simulation results show that, LEF cause significant changes in the operation states and the route-to-chaos of the lasers. We 
find that the values of the feedback rate at which the transition from CW to periodic oscillation or coherence collapsed or 
chaos state occurs decreases, as the LEF is increased. We identified the operation states of SL into six distinct operating 
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regimes, namely, CW, PO, PD, before route-to-chaos and 3xP and 4xP route-to-chaos and coherence collapse or chaos 
depending on the value of the LEF and OFB strength. The route-to-chaos is 4xP when FEF is less than or equal to three. 
The route is 3xP when FEF increases up to 7. The OFB strength when the laser transits from CW to PO or PD or chaos 
state decreases with the increase in the FEF. When the SFs operated under OFB and by decreasing the values of the FEF, 
the operations of the laser change from chaos operation to PO or CW operation, which is more stable and desirable. At 
higher values of FEF and when the injection current be near the threshold the laser stability is improved and moves toward 
stable operations. Finally we expect that SF subject to OFB exhibits much more stability under lower or higher values of 
FEF depending on the values of the injection current. 
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